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1. INTRODUCTION 


Because of the effects of material properties, structures 
constructed from orthotropic materials often require extensive analysis, 
even when they are geometrically simple. These analyses are commonly 
done using finite elements, boundary collocation, and Fourier analysis. 

The finite element method is well -documented in books such as 
Gallagher [1] and Zienkiewicz [2] and will be dealt with here only 
briefly. Finite element methods require the division of the region to 
be analyzed into smaller regions, or "elements." Within each of these 
regions, the form of an approximate solution for a stress field is 
assumed. Enforcement of boundary conditions and interelement continuity 
of displacement then gives a piecewise continuous solution for the 
stresses over the region. 

Similar to the finite element method is boundary collocation, or 
•point matching," which is described in the papers by Slot and Yalch [3] 
and Hooke [4]. For this method, an "exact" solution for the stresses in 
a region of materials is developed by requiring the stress field to 
match known values at given points on the boundaries. Unlike the 
finite element method, however, the stress field solution is good over 
the entire region being analyzed, not merely on one small portion of the 
region. 

Still another way to handle boundary value problems is to specify 
a solution form that allows the boundary conditions to be met exactly, 
rather than merely from point to point. This can be done with infinite 
series of polynomials, but is more often done with Fourier series. It 
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is to a form of Fourier series solution that this report is devoted. 

Timoshenko and Goodier [5] describe a fairly standard procedure for 
stress analysis using Fourier series. A Fourier series stress function 
is specified so as to match the conditions on two opposing boundaries 
and stresses are found from this function. A similar solution can be 
performed for the other two boundaries of the rectangular slab and, due 
to linearity, the stresses from the two solutions can be added together. 
A problem arises with this method, however. Since the two supposedly 
stress-free boundaries frequently acquire "residual" stresses from the 
solution, a new series solution must be formulated to cancel these 
stresses, i.e., to restore the stress-free condition on the adjacent 
boundaries. These series, in turn, will leave residuals on the original 
boundaries, and must in turn be cancelled. The process can be repeated 
as many times as necessary to reduce the residuals to acceptable magni- 
tudes. In 1944, Pickett [6] proposed a method to effectively "automate" 
this procedure. By taking a series expansion of the stress function in 
two dimensions, all boundary conditions can be met simultaneously. 
Details of the method used are described in Little [7]. Theocaris and 
Dafermos [8] have used Pickett's method to analyze a thin strip of 
isotropic material subjected to dislacement along a single boundary. 

In this report, a solution for the general plane stress problem in 
a rectangular slab of ortnotropic material is presented. The principal 
axes of the material are aligned parallel with the sides of the slab. 
Only stress boundary conditions are treated. However, the inclusion of 
displacement boundary conditions in the analysis is discussed briefly 
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in Appendix D. The solution Is valid for general boundary stresses, as 
long as they can be represented by Fourier series. The solution method 
is an extension and generalization of Pickett's method [6] tc allow the 
consideration of specia''ly orthotropic materials and general stress 
boundary conditions. The expansion of the capabilities of Pickett's 
method is made feasible by the existence of high speed computers which 
can quickly and easily handle large matrices and complicated algebraic 
expressions. 

When this solution is used, no redefinition of a finite element 
mesh is required to change loading conditions. Also (and probably of 
the most practical use in the long run), if regions over which solutions 
of this form are valid can be joined together to form more complicated 
shapes, they become "exact" elements in a new sort of finite element 
solution. Accurate results can then be obtained with far fewer elements. 

The solution generated is used to work two problems similar to 
tests used for material property evaluation. The stress fields found 
are examined for the effects of orthotropy, and possibilities for future 
use of the solution are discussed. 

1 .1 Method of Solution 

1.1.1 Development of Governing Equations 

The equation governing the plane stress distribution in an aniso- 
tropic slab is given by Lekhnitskii [9] as 


\ 
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^22 ^ - 2S«. 

3x"^3y 


+ (2S^2 ■*“ ^66^ 


_A 

3x^3y 


I - 2S^g 


3^(t) 

3x3y 


I 



( 1 ) 


where the stress function 41 is defined by the following relations: 




and T 


xy 


3^4> 

3x3y 


(2) 


The elastic compliances S. . are defined in Jones [10] and occur in the 

* J 

constitutive relations as 


S * ^ll'^x ^12‘^y ^16^ 

^y * ^12‘^x ^22^^y * he^xy 

Y = Si,-a + Sm-o + Sr/-T 

'xy 16 X 26 y 66 xy 

Positive values of the stresses a , a , and t are shown in Fig. 1. 

X y xy 

If the problem is specialized to descrit-: the stress distribution 
in a rectangular slab of orthotopic material where principal material 
axes coincide with the slab boundaries, the compliances S^g and $20 
vanish and the governing equation ( 1 ) reduces to 

( 4 ,) 

Using the standard engineering constants, this equation is 
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For the special case of isotropic material such that 


El - E2 and G12 - 2(1 +^Vi"2) ’ 

the governing equation reduces to the well-known bi harmonic equation for 
determining the stress function in plane stress: 


I? 


+ 2 


4 4 

a . a ct> _ 

2 4 ■ 

ax ay 3y 


0 , 


(5) 


The development of a solution to the orthotropic slab will be 
restricted to those problems where known stress distributions are applied 
on the boundaries x » ±£ and y = ±t, where i and t are the half-lengths 
of the slab in the x and y-directions, respectively (see Fig. 2). Thi*: 
restriction leads to the following boundary conditions: 


aJ^(^.y) = ri(y) 


aj^(-i.y) * r2(y) 

Txy(-^.y) = 

Oy(x,t) * r3(x) 

txy(x.t) = r^Cx) 

Oy(x,-t) * r^(x) 

Xjjy(x,-t) = rg(x) 


where r, are the given boundary stress distributions. These functions 
are not independent but from consideration of overall equilibrium of 
forces and moments are related as follows: 
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t 


I 


J [rT(y) - r 2 (y)] dy + 

-t 


[rs(x) - rg(x)] dx = 0 

-SL 


I t 

Lr 3 (x) - r^(x)] dx +[ [r^(y) - rg(y)] = 0 


-i 


(7) 


t 

[rgCy) + rg(y)]i dy 
-t 


i 

* 

[r^(x) + rg(x)Jt dx 
-I 


+ 


t 

[r 2 (y) - r^(y)]y dy + 
-t 



- r 4 (x)]x dx = 0 


Furthermore, the boundary stress distributions r^ must be expressible as 
Fourier series. 

Equation 4 and the boundary conditions of Eqn. (6) fully describe 
the governing set of equations to determine the distribution of stress 
in a rectangular slab of dimensions i, by t subject to any stress distri- 
bution on the boundary. 


2. SOLUTION TECHNIQUE 


A general solution of Eqn. (4) and (6) can be expressed as the 
tri Qonometri c expansion 

= c^x^ + C2xy + C3y^ + c^x^y + Cgxy^ 

dQ 

+ I cos(a y) + g (y) cos(B_x)] 

n»0 

00 

+ I [h^(x) sin(a^y) + k^(y) sin(S^x)], 
n=1 

where the polynomial terms are necessary in the solution procedure that 
follows to assure a general solution set, and where 


Ql 


n 



and 


mr 

I ' 


The functions ^p(x). 9n(y)> determined from the 

set of equations found upon substitution of the solution of Eqn. (8) 
into Eqn. (4). The result is a set of fourth-order ordinary differen- 
tial equations with constant coefficients in x for fp(x) and h^(x) and a 
corresponding set of equations in y for 9p(y) and k„(y)- These equa- 
tions admit of solutions of the form 




sx 


(9) 


which upon substitution in the ordinary differential equations all give 
the characteristic equation 


s^ - As^ + B = 0 , 


( 10 ) 
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where 


A “ — 



and 



(10) 


The roots to the characteristic equation are functions of material 
properties only. There are three possible solut'on forms; these are 

1) roots are real, in pairs — s = ±u, ±v (u,v > 0, u v) 

2) roots are real and equal in pairs s = ±u (u > 0) 

3) roots are complex — - s = ±u ±iv (u,v > 0; i = /-T) . 

These cases are now explored separately. 


CASE 1 

When the roots are real in pairs, case 1, the functions ^^(x), 
g„(y), >^n(y) ^^ave the form 


uoi_x va_x va X 

* W * ^31* * W 


u 8 _y v 8 y -ue y -vB y 

g„(y) “ a5„e * ^81.® 

ua X va X 

h„(x) ■ a,„e t a,„„e ♦ a,,„e * a,2„e 

u8 y v8„y -uB V -v8 y 

" 'I3n® * *14n' * *15n® * “l6n® 


( 11 ) 


To take advantage of symmetries which might occur in the formulation of 
specific problems, the solution form of Eqn. (11) is recast in terms of 
even and odd functions to yield 
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* [SSl^ cosh(ua^x) + ASl^ sinh(ua^x)]/cosh(ua^£) 

+ [SS3^ cosh(voi^x) + AS3^ sinh(va^x)]/cosh(va^t) , 

9n(y) = [SS2j^ cosh(uB^y) + AS2^ sinh(uB^y)]/cosh(u6^t) 

+ [SS4^ cosh(vB^y) + AS4^ sinh(vB^y)]/cosh(v8^t) , 

( 12 ) 

» [SAl^ cosh^ua^x) + AAl^ sinh(ua^x)]/cosh(ua^£) 

+ [SA3^ cosh(va^x) + AA3^ sinh(va^x)]/cosh(va^£) , 

kp(y) = [SA2^ cosh(uBj^y) + AA2^ sinh(uB^y)]/cosh(uB^t) 

+ [SA4^ cosh(vB^y) + AA4^ sinh(vBj^y)]/cosh(vB^t) . 


The terms cosh(uoi^O. cosh(va^O» cosh(uB^t), and cosh(vB^t) are 
inserted to avoid numerical instability in computation for large values 
of n, and the coefficients (SSl^, AS2^, etc.) are named in such a way as 
to denote states of symmetry in x and y. The first letter in eacn coeffi 
cient describes the symmetry with respect to the x-axis while the second 
letter denotes symmetry with respect to the y-axis ("S" denotes symmetry, 
while "A" denotes anti -symmetry) , Roots corresponding to case (1) occur 
in the analysis of highly unidirectional materials such as fiber- 
reinforced composites. In fact, any orthotropic material in which either 
of the normal moduli E-j and E 2 is significantly larger than the shear 
modulus G ^2 class. 


CASE (2) 

lihen the roots are real and equal in pairs, as occurs for isotropic 
materials, the functions fp(x)t 9p(y)* ^n(^)* similarly have 


the forms 


12 


^^(x) = [SSl^ cosh(ua^x) + SS3^ ua^x sinh(uoj^x) 

+ ASl^ sinh(ua^x) + AS3^ uo^x cosh(uo^x)]/cosh(uotj,Jl) t 

9n(y) * [SS2^ cosh(ue^y) + SS4^ u6„y sinh(u6^y) 

+ AS2^ sinh(u6„y) + AS4^ u6„y cosh(u6j,y)]/cosh(uBpt) 

(13) 

h^(x) » [SAl^ cosh(ua^x) + SA3j^ ua^x sinh(ua^x) 

+ AAl^ sinh(ua^x) + AA3^ ua^x cosh(ua^x)]/cosh(ua^£) 

kj,(y) = [SA2^ cosh(uB^y) + SA4^ uS^y sinh(uS^y) 

+ AA2^ sinh(uS^) + AA4^ uB^y cosh(uBj^y)]/cosh(uBj^t) . 

As with the previous case, the coefficients denote synmetry conditions. 

CASE 3 

When the roots are complex, case (3), the functions ^^(x), 9p(y)> 
h^(x), and k^(y) have the forms 

f„(x) = [(SSl^ cos(va^x) + ASlj. sin(va^x)) cosh(ua^x) 

+ (SS3^ cos(va^y) + AS3^ sin(va^x)) sinh(ua^x)]/cosh(uaj^£) , 

9n(y) = [(SS2^ cos(vB^y) + AS2^ sin(vB^y)) cosh(uB^y) 

+ (SS4^ cos(vB^y) + AS4^ sin(vBj,y) )sinh(uB^y)]/cosh(uB„t) , 

(14) 

h^(x) » C(SA1^ cos(vaj^x) + AAl^^ sin(va^x))cosh(ua^x) 

+ (SA3^ cos(va^x) + AA3^ sin(va^x)) sinh(ua^x)]/cosh(uo^£) 

k^(y) '= C(SA2^ cos(vB^y) + AA2^ sin(vB^y)) cosh(uB^y) 

+ (SA4^ cos(vB^y) + AA4^ sin(vBj^y)) sinh(uBj,y)]/cosh(u6^t) . 


13 


As with the previous two cases, the coefficients denote symmetry condi- 
tions. While case (3) is not fully developed in this report, it is 
nevertheless quite significant. Any material with large shear modulus 
G.J 2 (as compared to the normal moduli and E 2 ) fits into this case. 

For instance, graphite/ epoxy in unidirectional form has material proper- 
ties [n] 

E^ = 19.2 X 10® pxi 
E 2 “ 1.56 X 10® psi 
G .,2 * 0.82 X 10® psi 
V-J2 * 0.238 

and falls, into case (1). However, using this material in angle-ply 
laminates having the stacking sequence C-^lsymmetric ^*^® shear 

modulus while lowering the maximum normal modulus, and gives properties 
corresponding to case (3) when 0 lies in the range 

25. r < 9 < 64.9°. 

Since laminates of this form are used frequently in composite structures, 
this latter case obviously is important. The purpose of this paper is 
to describe a general technique for solution and the general development 
of equations will be limited to cases (1) and (2). The procedure for 
incorporation of material systems requiring case (3) solutions follows 
directly using the same development as described for the first two cases. 




h 
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Reformulation of the Stress Function 

The equation for the stress function <(), Eqn. (8), can be recast to 
take advantage of even and odd functions as follows: 

^ “ '^SS ^SA ^ ^AA 

where 


"SS 


jR 


xss^ * 2 ^yss^ 


I “*<V> ^ I 


n=0 


n=0 


♦SA • J * I ’‘nsa * I 


n»l 


n=0 


00 00 

♦as ' I "xas*)' * nio *nas '“'V’ * I Vas 


n=l 


♦aa ' "aa*^ * I ’'naa ''"'V) * I ''naa 
n*l n=l 


In these expressions, the terms and X^^^^ are functions 

of X, while the terms Y^^^, and Y^^^ are functions of y. 

The forms of these functions depend on the roots to the governing 
differential Eqn. (4), as discussed previously, and are presented later. 
The terms R^^^, R^^^, and R^^ are arbitrary constants. 

Using the definition of the stress function (Eqn. (2)), stresses 
are derived from the stress function by differentiation. These stresses 
are given in terms of the stress function in Appendix A. 
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For Case (1), roots real and in pairs, the functions etc. 

are given by 

» SSl^cosh(uaj^x)/cosh(ua^jl) + SS3^cosh(va^x)/cosh(va^t) 

Xnsa ’ SAlj^cosh(ua^x)/cosh(uaj,a) + SA3^cosh(va„x)/cosh(va^l) 

^nas * ASlj^s1nh(ua^x)/cosh(uaj,Jl) + AS3j^sinh(vOj^x)/cosh(va^t) 

Xnaa “ AAl^s1nh(ua^x)/cosh(uo^l) AA3y^s1nh(vo^x)/cosh(va^t) 

^nss ” SS2^cosh(u6j,y)/cosh(u6j,t) + SS4^cosh(vB^y)/cosh(vS^t) 

Ynsa “ SA2j^s1nh(uB^)/cosh(u6j,t) SA4^sinh(vg^y)/cosh(vBj^t) 

■ AS2^cosh(u6jP')/cosh(u6^t) + AS4^cosh(v6^y)/cosh(v6pt) 

Y„aa = AA2„s1nh(u6„y)/cosh(u6„t) + AA4„sinh(vB v)/cosh(v8_t) . 
ndo n n n n n n 

For case (2), roots real and equal in pairs, these functions are given by 

^nss " ^^^^n ^ sinh(ua^x)]/cosh(ua^S,) 

X^gg = cosh(ua^x) + SA3^uo^x sinh(ua^x)]/cosh(ua^£) 

X^^^ = cosh(ua^x)]/cosh(ua^t) 

^naa * ^^^n siih(uoi^x) + AA3^ua^^x cosh(ua^x)]/cosh(uot^il) 

Ynss = [SS2j^ cosh(u8^y) + SS4^uB^y sinh(uSj^y)]/cosh(u6j^t) 

Y^sa “ tSA2^ cosh(u6„y) + AS4^u8^y sinh(u8^y)]/cosh(u6^t) 

Y«ac * CAS2„cosh(u6 v) + AS4 u8_ysinh(u6 y)]/cosh{u8_t) 
nds n n n n n n 

^naa * ^^^n "*■ ^^n“®n^ cosh(u6^y)]/cosh(u8j,t) . 

Finally, for case (3), complex roots, these functions are given by 

^nss ' [SSl^ccs(vd^x)cosh(uaj^x) + SS3^sin{vapX)sinh(uapX)]/cosh{uaj,e) 
Xnsa ’ [SAl^cos(va^x)cosh(uaj^x) + SA3^sin(vo^x)sinh(ua„x)]/cosh(uaj^i) 
^nas “ [ASlpSin(vapX)cosh(uOj^x) + AS3^cos(va^x)sinh(uaj^x)]/cosh(ua^O 
Xnaa “ CAAlnSin(vanX)cosh(uanX) + AA3nCOs(vanX)sinh(uanX)]/cosh(uoin)l) 


1 
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Yfiss “ CSS2^cos(v6j^y)cosh{u6j^y) + SS4^sin(v6j,y)slnh(u6j,y)]/cosh(u0j^t) 
Y^sa “ [SA2^sin(v6^y)cosh(u6^y) + SA4^cos(v8j,y)sinh(u6^)]/cosh(u8^t) 
^nas * [AS2^cos(v8^y)cosh(u8j^y) + AS4^sin(v8^y)sinh(u8^y)]/cosh(u8„t) 
^naa * CAA2^sin(v8^y)cosh(u6j,y) + AA4^cos(v6^y)sinh(u8^y)]/cosh(u8^t) 

08 ) 

Formulation of Boundary Conditions 

The stress function 41 as expressed by Eqn. (15)- (18) is composed of 

three types of functions which contain unknown coefficients. These are 
1 2 

the polynomials ( 2 " the even series functions ^nsa’ 

''nss' ''nas>- 

Each of the series functions contains two unknown coefficients, so the 

even series functions contain eiaht (SSI , SS3„, , AS2„, AS4„) for 

each value of n, and tha odd series functions contain eight (ASl^, AS3^, 

, AA2^, AA4^) for each value of n. Thus, 16 unknown coefficients 

must be found for each value of n in the series, and another five must 

be found to account for the polynomial terms (R„,,. R„c,» 

xss y^s yso xas 

"aal- 

To obtain the solution for the coefficients in Eqns. (15)-(18), the 
boundary conditions (Eqn. (6)) are expanded in general Fourier series, 
with the form 

go 

^/(i.y) * 2T ^ 1 q ^ T ^^1 h 

t 

Cl = f r,(y) dy 

0 it 


where 
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n 


t 

» 

r^(y) cos(a^y) dy 
-t 



t 

' 

r^(y) sin(a^y) dy . 

/ 

-t 


Similar expansions hold for the other applied boundary stresses. See 
Appendix B for these expansions and their Fourier coefficients. 

To take advantage of possible symmetries which might occur with the 
stresses applied on the boundaries, boundary conditions are obtained by 
equating the sums and differences of the applied stresses (given in 
Appendix B) at the boundaries to the corresponding sums and differences 
of the general expressions for stress (see Appendix A) evaluated at the 
boundaries. The resulting equations are then manipulated to extract 
expressions for the arbitrary coefficients in the following fashion. 

For example, on the boundaries x = ±£, 

r^(y) - T2(y) = Oj^(£,y) - o^(-£,y) (20) 

Substituting for r^(y) and '^' 2 ^^) expanded in Fourier series (as given in 
Appendix B) and for as given in Appendix A and evaluated at the 
boundaries, this equation becomes 


2t ^l-2(, t Jj^l-2„ ^l-2„ sin(a^y)] 

u n* 1 n n 


^'^xas ■ ^ ^Iq °‘n^nas^*'^ cos(ci|^y) 




( 21 ) 



18 


where terms of the form C^_2 are a shorthand notation for (C-j -C2 ). 

P P P 

Recall that terms of the form are functions involving unknown 

coefficients such as ASl^, evaluated at 2 .. To relate the unknown 

coefficients to the known Fourier coefficients of the applied stresses, 

both sides of the equation are multiplied by either sin(a^y) or cos(a^y) 

and integrated along the boundary from -t to t (for stresses on the 

adjacent boundaries, the boundary condition equation is multiplied by 

sin(B^x) or cos(S^x) and integrated from -2 to 2 ). Continuing this 

process for the remaining boundary condition relations, 

f'^(y) + r2(y) » + a^(-2,y) 

rjCx) + r^(x) = ay(x,t) - 
r3(x) + r^(x) » Oy(x,t) + a^(x,-t) 
rs(y) - rg(y) • T^^(Z.y) - T^/-t.y) 

T5(y) + rg(y) = T^y(^.y) + 

r7(x) - rg(x) = T^y(x,t) - T^y(x.-t) 

r7(x) + rg(x) = ’ 

the following set of equations result for the unknown coefficients with 
all material cases: 


SS 1 p * (ElAp) i (ElBp„)SS 2 „ = (ElAp){(ElCp) + j (ElDp„)) 
n=l n=l 


(FlAp) I (F’Bp„)SSl„ * SS 2 p = (nApH(FlCp) * I (FlDp„)) 

n=l n*l 
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SAlp ^ (E2Ap) I (E2Bp„)SA2„ - (E2Ap)( £2Cp) > ^ (E2DpJ} 
n»l n*1 


(f2Ap) I (F2Bp„)SAl„ + SA2p • (F2Ap){(F2Cp) + I (F2Dp„)} 
n=l n=l 


W’p + (E3Ap) I (E3Bp„)AS2„ ■ (E3Ap){(E3Cp) t J (£30^^)) 


( 23 ) 


n»l 


n*l 


(F3Ap) I (F3Bp„)ASl„ + AS2p = (F3Ap)((F3Cp) + I (F3Dp^)} 


n»l 


n=1 


AAlp t (E4Ap) I (E4Bpp)AA2„ = (EA4p)t(E4Cp) » I (E4Dp„)) 


n=l 


n*l 


(F4Ap) I (F4Bp^)AAl^ + AA2p = (F4Ap){(F4Cp) + I 0^4Dp^)> • 


n=1 


n=l 


Furthermore, equations of th§ following type also result (these examples 
are for Case 1 material): 


SS3„ = 


1 

'%-6p 

V tanh(vcipil) 

2 

L2«pt 

1 

^7-8p 

V tanh(vBpt) 



- SSI u tanh(ua ?,) 


SA3_ = 


-1 


p V tanh(vap2,) 


+ SAl.u tanh(ua„J,) 

2a^t P P 

p 


SA4. 


tanh(uB„t) C-, , 

^^^p tanh(v6„t) 7^2 


p-' 26pJl tanh(vBpt) 


(24) 
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" tanh(vapA) 



AS4 


- 1 


rc 


p V tanh(va_t) 


^ + AS2 u tanh(u3 t) 




AA3. » 


- 1 


p tanh(vap£.) 


5 

AAl tanh(ua„£) + 

" ^ 2a?t 


“ tanh(vSpt) 


AA2^ tanh(u6„t) + 

P P 23^i 


For a complete definition of the terms in Eqn. (23), as well as for the 
Case 2 expressions corresponding to Eqn. (24), see Appendix C. 

The Eqns. (23) can be paired in matrix form according to symmetry 
conditions as, for example. 


r I 1 El 1 
1 

. 1 

mj 

FI T 1 

C 

1 

1 CVJ 
f to 
1 to 
1 

1 CD 
1 r— 
1 O- 


n 


n=l,2, — (25a) 


where, after partitioning, [I] is the identity matrix, [El] and [FI] are 
square matrices dependent on material properties and slab geometry; 
{PIA}^ and [P1B}|^ are column vectors dependent on material properties, 
slab geometry, and boundary conditions; and (SSl}^ and {SS2}^ are column 
vectors containing the unknown constant coefficients of the series 
stress function. Similarl>, 


E2 


SAl 


F2 


! I 


SA2 



n » 1,2, — (25b) 
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I 1 
1 
1 

E3 

( ASl 


P3A 

F3 i 

I , 

( AS2 

1 

P3B 

I i 
1 
1 

E4 ■ 

j AAl 


P4A 

F4 i 

I 

( AA2 

n 

P4B 


n * 1,2 -- (25c) 


n = 1.2, — (25d) 


Each set of Eqn. (24) may be solved Independently; for example, solving 

for (SS1}„ and {SS2}„ gives 
n n 

SS2 ^ = [[I] - [F1][E1]] {PIB} - [F1]{P1A}^ 

SSI ^ = PIA - [E1]{SS2}^ . 


Using the unknown coefficients determined from the above matrix solutions, 
the remaining unknown coefficients can be found from Eqn. (24) for Case 
1 or from the corresponding equations for Case 2 material. 

The only terms which are still unknown at this point are the 
coefficients of the polynomial terms in the expression for the stress 
function (Eqn. (15)), These are found in similar fashion to the above 
using the boundary relationships of Eqns. (20) and (22) and n = 0. 

The polynomial terms R^,., and R^._ are then given by 

ySo XaS 


'3-4, 


ysa 


xas 


4U 


( 27 ) 


i 


i 

I 



for both materials (1) and (2). The polynomials 

xss yso QQ 

dependent^ not only on material properties, but also on the solutions for 
the coefficients SSl^^, SS2^, etc. For material case 1, roots real and 
in pairs, these terms are given by 

*^xss “ ^ ^1+2 ‘ I r tanh(u6„t) SS4^v tanh(v6„t)] 

° n-1 

%ss * ^ ^3+4 ■ ^ tanh(ua^il) + SS3j^v tanh(vaj^t)] 

° n*l 

^5+6f. « 6 

*^aa ^ " “?t ^ ~T ^^n tanh(v3^t)] . 

n*l 

For material Case 2, roots real and equal in pairs, these polynomials are 
given by 

1 “> u8 

\ss=W'u2 ■ IH)" [SS2„ tanh(u6„t) 

^ n=l 

+SS4„{tanh(u6„t) + u6„t}] 
n n n 

- I H)"t! 1[SS1„ Unh(ua^) 

“ n.l 

+ SS3„{tanh(ua t) + ua„t}] 
n ' n n ■' 

^5+6n « S„ 

«aa = - ^ ^ T tanh(uB„t) . AA4^uB„t] . 

n»l 

Using Eqns. (23) th»*ough ) and Appendix C all of the unknown 
coefficients can be solved. Then, using the equations given in Appendix 
A, the stresses at any point in the material can be found. 



3. COMPUTER IMPLEMENTATION 


To solve the equations described previously, a computer program was 
written and programmed on a Control Data Corporation Cyber 173 computer 
at the Langley Research Center of the National Aeronautics and Space 
Administration. A brief description of this program follows. 

After receiving input of material properties and the Fourier series 
coefficients for the boundary stresses, the program calculates the roots 
of the characteristic equation from the material properties and decides 
which solution form to use. The subroutine corresponding to that form 
is then called to solve for the unknown coefficients. The various 
symmetry cases are solved one at a time by computing the elements of 
matrices E and F as required with the values indicated in Appendix C. 

A matrix manipulation subroutine is then called to perform a set of 
operations common to all symmetry cases. At the end of these operations, 
the coefficients for the symmetry cases in question are put into the 
proper storage locations and the matrix routine returns control to the 
routine which fills the matrices. The next symmetry case is then solved 
in the same fashion, until all symmetry cases are done. Any single type 
of symmetry can be considered (to the exclusion of all others), or the 
program can solve all four types in turn, thus solving a general 
(asymetric) loading case. A check on the accuracy of the solution is 
given by printing out the input boundary conditions (found by summing 
the input Fourier series) in columns beside the corresponding output 
boundary conditions (found by summing the "exact" solution on the 
boundaries). This feature can also be used to check for errors in input. 


'£»■ 


23 
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Finally, the locations for which stress values are desired are checked, 
stresses are calculated at these locations, and a subroutine Is called 
to print or plot these stresses in useful form. 

Several problems arise in implementing the computer code. It has 
already been pointed out that terms based on the maximum values of x and 
y and on the value of n are included in the solution. These terms are 
included to prevent terms involving hyperbolic functions from becoming 
too large and destroying the accuracy of the series (by requiring the 
differences of large numbers to be much smaller than the numbers them- 
selves). For plates with high or low aspect ratios (S,/t greater than 
about 4 or less than about 0.25), the arguments of terms such as 
cosh (uot^Z) become too large for the computer to handle when n becomes 
large (n approximately equal to 50). Thus, the solutions for these 
plates require that fewer terms of the series solution be taken. It is 
possible that this may lead to unacceptable accuracy if the loading is 
such as to require a large number of terms for convergence of its 
Fourier series representation. 

Because of the nature of Fourier series expansions, convergence 
problems might occur for some types of loadings at the corners of the 
slab. 

3.1 Sample Problem Solutions 

Two practical sample cases were chosen as examples of potential 
applications of this solution procedure. The sample cases presented 
were selected for their similarity to actual problems encountered in 
materials testing and to demonstrate the wide applicability of the 



25 


program developed. 

Solutions were generated using 35 and 50 terms in the series formu- 
lations. Plots were generated taking cuts at constant x/e, equal to 0.0, 
0.2, 0.4, 0.6, and 0.8 with -t < y < t. In all cases, t * 1.0, while 
£ » 1.0, 2.0, 4.0 and 6.0. Thus the slabs analyzed had aspect ratios 
A.R. “ l/t ■ 1.0, 2.0, 4.0 and 6.0. 

The analysis also involved two different materials. One of these 
was isotropic with the properties of steel; 

= E 2 = 30.0 X 10® psi 

G .,2 = 12.0 X 10® psi 

v-|2 = 0.25 . 

The other material was orthotropic, having the properties of uni- 
directional graphite/epoxy from [11]: 

E^ = 10.2 X 10® psi 

E 2 = 1.56 X 10® psi 
G^2 = 0-82 X 10® psi 
v^2 “ 0.238 . 

The loadings, to be referred to subsequently as the "tensile coupon" 
and "rail shear" cases, are described as follows. 
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! 

I , 


3.1.1 "Tensile Coupon: Load Case 

The "tensile coupon" test has the following boundary conditions: 

T (x,±t) » 0 for -2. < X < -(0.95)il 
xy • 

-(0.90)Jl < X < (0.09)4 
(0.95)4 < X < 4 

= ±(20)t/4 for -(0.C£‘, !. < X < -(0.90)4 

= ±(20)t/4 for (0.90)4 < X < (0.95)4 

= 0 

c^j((ii.y) = 0 

(7y(X,±t) = 0 . 

These boundary conditions (see Fig. 3) are chosen to give a very crude 
approximation of the stresses on a rectangular coupon clamped in a test 
rig and being pulled for the determination of tensile strength, ihe non- 
zero values of boundary shear are selected to give an average normal 

stress C7^ of magnitude 1.0 when the applied stresses have diffused 
completely. 

3.1.2 "Rail Shear" Load Case 

The "rail shear" test has the boundary conditions 

T (x,+t) * 1.0 for -(0.80)4 < X < (0.80)4 
xy - 

= 0 for -4 < X < -(0.80)4, (0.80)4 < x < 4 

a„(x,±t) = — yx for -(0.80)4 < x < (0.80)4 

y (0.64)4^^ ■ ■ 

= 0 for -4 < X < -(0.80)4, (0.80)4 < x < 4 

T^y(±4,y) =• 0 
aj^(±4,y) = 0 . 


I 
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Plot of T for Isotropic " 
xy 


Fig. 14 


Rail Shear Specimen" of A.R. = 6 
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"rail shear" specimen. All specimens showed essentially identical 
behavior, with only the magnitude of the stress changing. This is an 
aspect ratio effect, since higher magnitudes of o^ are required to 
satisfy moment equilibrium for lower aspect ratio specimens. 

Aspect ratio provides most of the variation in normal stress 
distribution for the "rail shear" tests. The distributions of in 
orthotropic and isotropic plates have nearly identical forms and magni- 
tudes for a given aspect ratio, though the maximum magnitude of is 
slightly lower for orthotropic materials at low aspect ratios. At 
high aspect ratios, the difference in magnitudes due to material 
properties is negligible. See Figs. 11 and 12 for plots for an 
isotropic plate of aspect ratio of 1.0 and an orthotropic plate of aspect 
ratio 6.0, respectively. 

Figures 13 and 14 show the shear distribution in isotropic plates of 

aspect ratios 1.0 and 6.0, respectively. A comparison of these plots 

serves to illustrate tne smoothing of the curves for x as aspect ratio 

xy 

increases in both isotropic and orthotropic plates. Orthotrcpic plates 

of low aspect ratio, however, tended to display i behavior similar 

xy 

to isotropic plates of slightly higher aspect ratio (see Fig. 15). The 

material properties were of little consequence to the shear stress field 

at higher aspect ratios. At the edge of the stress boundary {x/l * 0.8), 

the behavior of y = ±t shows a value of i = 0.5 instead of the 

xy 

anticipated value of 1.0. This is a result of the fact that Fourier 

series give average values at discontinuities, so the applied boundary 

stress input at this point is, in fact, x = 0.5, rather than 1.0. 

xy 
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3.3 Conclusions and Recoirwiendations 

A method for solving exactly the stress field in a slab of ortho- 
tropic material was developed. Tests were run to determine the useful- 
ness of the method. The procedure suggested by Pickett [6'j shows consi- 
derable promise for use in the analysis of both isotropic and ortho- 
tropic materials. Accurate results were obtained with fairly low run 
times. 

One of the expected benefits of the method proved to be nonexistent. 

It was originally anticipated that all rectangular membranes could be 
analyzed with similar accuracy, but this was found not to be true. The 
inability of the computer to handle large arguments of hyperbolic sine 
and cosine functions limits the number of terms which can be taken in 
the series solution. For plates with very high aspect ratios, this 
limitation can become critically small and destroy the accuracy of the 
solution. 

The first and most important task for future researchers using 
Pickett's method is to develop the solution for Case 3, involving complex 
eigenvalues for the governing differential equation (Eq. (4)). Due to 
the significance of materials which fall into this case (for example a 
^"^^°^symmetric graphite/epoxy) the inclusion of this 

solution is important. Once this is done, a thorough study of convergence, 
accuracy, and computing cost for Pickett's method needs to be undertaken. 
These parameters should be compared to similar parameters for one or 
more finite element programs. 
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Assuming that the suggested research proceeds with no major 
difficulties, a displacement solution should be formulated following 
the approach suggested in Appendix 0. This would allow the solution of 
more general problems by allowing completely general boundary formu- 
lations. 

Once a displacement solution has been developed, a method for 
joining one of these rectangular slabs to another needs to be developed. 
The slab for which the displacement solution holds could then be turned 
into a finite element, capable of deforming compatibly with its neighbors 
at every point of their common boundaries. For the time being, this 
joining of solutions must necessarily involve a minimal number of 
elements, due to computer limitations. However, more powerful computers 
already being developed make ihe concept of "exact" finite elements 
appear feasible in the very near future, and the technology required 
for them needs to b-; developed now. 

An attractive feature of using Pickett's method for finite elements 
is that only a few elements may be required even when high local stess 
gradients exist. Presently, large numbers of elements are required for 
accurate solutions in such regions. 


k 
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APPENDIX A 


STRESSES IN TERMS OF THE STRESS FUNCTION 


Using Equation 15 and the definition of the stress function (Equation 
2), the following relations for stresses can be defined: 


a * a +a +o +o 
X xss xsa xas xaa 


0*0 +0 +.7 +0 

y yss y^a yas yaa 


T=T +T +T +T . 

xy xyss xysa xyas xyaa 


(A-1) 


where 


0 * — ^ 6 
xss .. 2 ^ss 


r 


xyas 3x3y ^as 


» etc. 


With these ter"S defined, the following component stresses for 
result: 


0 = R - ? ot^X cos(cx y) + T Y cos(B x) 

xss xss ^“Q^*n nsj ^ n*0 ^ 


* ’'xas* - „;„=‘^nas"’®<V> * 


(A-2) 


- ? 


- ^n^naaS''n(a„y) + ® Y„^,Sin(3„x) 

xaa ^ I n naa n , naa n 

n= 1 n=i 


i 
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Similarly, for 0 ^: 

”ysa = ''ysa^ * 

Vs ° „Vms“*<V> - J e^nas"'"<V> 

n*u n- I 

"yaa ' „|,Ca'"<V) ' J/^naa**"*^*) ' 

Finally, for t : 

^yss ’ J,”'n’'nss^*"<V' * J/n^ss*'"* V> 

Vsa = -"ysa’' ' J/a*nsa“*<V> * J,8n’'nsa^''"(«n*> 

Vas “ -"xasJ' * J,=‘n*nas='"<V> ' J/n’'Ms'='>*‘8„x) 

"xyaa ' -"aa ‘ ’ J,6n''naa">5(Sn">- 

In all of the above expressions. 


(A-3) 


(A-4) 


=dl' 


" d 
X = -^ X 

dx"^ 


Y = ~y 
dy 


dy^ 

The X and Y Terms are dependent on material properties and are defined 


in Equations 16-18. 
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APPENDIX B 

FOURIER EXPANSIONS OF BOUNDARY CONDITIONS 

The boundary conditions of Equation 6 can be expanded in general 
Fourier series. These expansions are given below. 


<J ,(t.y) • ItC, * ? i c,tos(v) ♦ S, Sln(v) 

n» I rj^ n 


Cl = i ^{y)dy 

I o * 


= / r^(y)cos(a^y)dy 
n -t 


= / r^(y)sin(a^y)dy 
n -t 


^ “ It 4o t 


C2 cos(a^y) + S^ sin(o^) 
n n 


^2o ^ 


^2 ^ / r2(y)cos(a^)dy 
^n -t 


53=/ r 2 (y)sin(a^y)dy 
n -t 


a y(x.t) 


2^ So ^ 


cos(S^x) + S^ sin(B^x) 


j r.,(x)dx 
-i ^ 
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C 3 » / r 3 (x)cos( 0 ^x)dx 


'n -i 


S3 = / r3(x)sin(6^x)(lx 


'n -I 


• k ^4, * 1 lV“‘V ^ 


C4 - j 'f.(x)dx 

a 

C4 = r^(x)cos(B^x)dx 


'n -il 


S 4 * J r^U)sin( 6 „x)dx 




\yU.y) • k ^5, ^ T 1=5 «=(V> ^ 


Cc “ / rcly)dy 


-t 
t 


Cc = / r5(y)cos(a^y)dy 


-t 


Sc = / rc(y)sin(a^y)dy 

-t 


T^y(-t.y) 


! 

■kk + ? I Cj cos(a y) Sg sin(a^) 

^t| 60 n=l ^ '’n n 

t 

^ Tg (y)dy 
t 

^ rg(y)cos(aj^)dy 


Sg = / rg(y)sin(a^y)dy 


'n -t 
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= k‘^7. * J, I [C7„«*<V> " 


= / r7(x)dx 


-I 


C7 = / r,(x)cos(8 x)dx 
n -l 


S7 » / r,(x)sin(S^x)dx 
'n -H 


Tjy(X.-t) = ^ Cg^ ■* J 


Cg cos(6„x) + Sg sin(3^x) 
°n n 


Cg^ = --g(x)dx 


Cg = / rg(x)cos(6^x)dx 
n -il 


Sg = / rg(x)sin(6^x)dx 


Due to linearity, expressions of the form 




1 

It 


(C, + C. ) + ? r 

lo t 


(C, + Cp )cos(a y) + (S, + $2 )sin(aj^y) 

n n n n 


and 




It ■ ‘'2o^ t 


(C^ - C2 )cos(a^y) + (S^ - $2 )sin(a^y) 

n n n n 
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are valid. These forms are also useful in the solution technique of this 
paper, so the following shorthand expressions are introduced: 


'' 1+2 ■ ^1 ^2 
n n ^n 


^1-2 “ ^1 * 4 

n n n 


^l+2„ “ ^1„ \ 

n n n 


^l-2„ “ ^1 ■ ^2 

n n n 


^ 3+4 ^ ^3 ^4 

n n n 


^3-4 " ^3 " ^4 
n n n 


^3+4 “ ^3 ^4 

n n n 


^3-4 ^ h ■ ^4 
n n n 


^5+6„ “ ^6 

n n n 


^5-6 ' * ^6 

n n n 


^5+6 ’ 

n n n 


^ 5-6 * ^5 ‘ ^6 
n n n 


^7+8 " ^ \ 

n n n 


^ 7-8 " ^7 ■ ^8 
n n n 


^7+8 " ^7 ^8 

n n n 


^7-8 " ^7 ■ ^8 
n n n 


Using these. 




1 r . s 1 

Jt ‘^1+2. * t 


C,*2 cos(s„j) + S,^2 sin(s„y) 


1 
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za 


‘^7-80 J, 


and so forth. 
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APPENDIX C 

EXPRESSIONS FOR SOLUTION TERMS 


The matrix expressions given by Equations 23 and 25, along with 
the explicit expressions for solution terms as in Equation 24, are the 
backbone of the solution technique presented in this paper. Each of 
the three material cases provides a distinct set of terms for Equation 
23 as well as distinct sets of equations with the forms of Equations 
24, 27, 28 and 29. The terms and equations corresponding to material 
cases (1) and (2) are given below: 

Case 1: Roots Real and in Pairs (±u, ±v; u,v>0; u ^ v) 

Symmetric x, symmetric y (SS): 

SSlp t (ElAp)J^(E16p„)SS2„ " (ElAp){(ElCp) * (EID^J) 




SS2p = (FlAp){(FlCp) 


n=l 




where 


E1A„ = 

P 2 


1 




1 - 


utanh(ua_£) 


vtanh(vapA) 




uv 


2 






tanh(uB^t) 


E'S = ■ 


^5-6p ^ c 

vtanh(vapt) ^l+2p 
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EID 


pn 


2(-i)"^Pa 


" i[(v 6 „£)^" + cxb 


FI A. 


1 


28^1 


1 - 


utanh(u6pt) 

vtanh(u6"t) 


F 18 p„ . - 4 (-l)"*Po^ 


u’ 

T~^- 


uv 

v 2 . rt 2 - 




FTC, 


- 7-8p ^ c 

vtanh(v$pt) ^3+4p 




Then, using these coefficients. 


^^^p “ vtanh(v« 


rs 


p"^ L 2 cijt 


5 ^ 

- SSI utanh(uoi £) 


^^^p ~ vtanh(v6„t) 


p'-' L 26 "i 

K 


- SS2_utanh(u6_t) 
r.C„ P D 


"xss ' 3t ‘^U2, • J, t^SS2„utanh(u8„t) ♦ SS4 


Vs ‘ k - V'”" rt'S1„uUnh(ua„l) » SS3„ 

n= I 


Symnetric x, antisymmetric y (SA): 


SAT + (E2AJ ? (E2B„JSA2„ = (E2AJ {(E2Cj + ? 


n*l 


pn' n 


n=l 


tanh(uaj^il) 


^vtanh(v6„t)] 
vtanh(va^i)] . 


1 
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where 


E2A 


1 


P 2 


1 - 


utanh(ua_Jl) 


vtanh(vapil)j 


E2% • 


( U 6 „)‘ 


( V 8 „)' 


C(u6„)^ + C.2] [(V6„)^ ♦ CL^: 


E2C = « 5 

vtanh(vapS) l+2p 


E20„, ’ 4(-l)P% 


(^8„) =3- 


4n 


'P" ■' ■' 'P C(v0^)2 + a^] Z&ll 


- 2 

tanh(u8pt) 

2Spt 

tanh(v8pt) 


F2Bp^ = 4(-l)P''%a2Sptanh(ua„i) 


1 


1 


C(%)^ " ^ 


^ ^7+8p ■*■ tanh(v8pt) ^3-4p ' gpt 


3-4o 


F2D„„ = 2(-1)P'^" K Cc 

P" t[(va^)2 + S^j 


S. 


Further SA relations are 


-1 


" vtanh(va_i) 


+ SAKutanh(ua_]l) 


L2a^t 


,!, ‘™pp” 


tanh(u6j,t) 


56 




’tanh(u 6 pt) 

tanh(v 6 pt) 


+ 


^ 3 - 4 p 

26^Uanh(v3pt). 


R = 
ysa 4 it 


Antisymmetric x. Symmetric y (AS): 

ASI + (E 3 A ) ? (E 3 B )AS 2 » (E 3 A ) {(E 3 C ) + ? (E 3 D„J) 
P P pal pn n p p pn 

(F 3 A ) ? (F 3 B JASI + AS 2 = (F 3 A ) {(F 3 Cj + ? (F 3 D_J} 

P P_] pn n p p p pn 


where 



1 


« 2 

tanh(uapS,)' 

Za t 

p 

^ tanh(vapi) 


E 3 B 


pn 


4 (-l)P'^"u 6 ^a tanh(u 6 t) [ ^2 T 


1 


. [(vB + ah 
p-* *^n' p-* 


E3C = S + ^ r . (-l)P 

^5+6p tanh(vcXpS,) h-2p ^ cxpi 


E 3 D„„ = 2 (-l)P'"" 5 5 - v -7 On 


F 3 A. = 


1 


26 ;^ 


utanh(uBnt) 
1 - P 

vtanh(v 0 _t) 


F 3 Bp„ = 




LC(%)^ ♦ Bp] 




tanh(ua i) 
n 
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Fir = ^ 7 - 8 p _ > 

vtanh(vBpt) ^3+4p 




Also, 


AS3_ - 


-] 

tanh ( vapil ) 


AST tanh{ua 5.) + 

^ ^ 2 a : 


1^1 

i»:t J 


AS4 


■1 


p " vtanh(v6_t) 


p ^' L 2 Sp)l 


+ AS2 utanh(u6.t)| 


5 , n-2o 

^xas 4 U ‘ 


Antisymmetric x, antisymmetric y (AA): 


AAlp t (E4Ap) (E4Bp„)AA2„ = (BAA^) {(£40^) + (E4D^„)) 


(F4Ap) (F4Bp„)AAl„ + AAB^ • (F4Ap) {(F4Cp) + {FAD^J} 


where 


E4A. 


2 “pt 


U - V 


tanh(ua„£) 

_g_ 


tanh ( vapi ) 


E4Bp„ » 4(-1)P^6, 


3 r 

" L [( u 6„)2 * 4l 


v 2 
,2 , 2 - 


[( V 6 „)‘ F apj 


tanh(uB^t) 


^^^p ~ tanh(vapJl) ^l-2p " ^5+6p 


k 
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(vS ) 

E4D = ^ ^ S- - 

"" * “p] 


F4A^ = — 

26^1 


1 


tanh(u8 t) 
‘ tanh(v6„t) 


F4B„ » 4(-l)P‘^"aJ 
pn n 


u2 




LC(ua„)^ Bp [(va„) .p 




tanh(ua^il) 


‘ tanh(vSpt) ^3-4p ” ^7+8p 


= 2(-l)'’'^ T- Pn 


Also, 


-1 


^^p “ tanh(va_il) 


AAl_tanh(ua_Jl) + — -y ^P - 
P P 0-2^ 




^^p “ tanh(vBpt) 


AA2_tanh(u6.t) + ^ - I^P 

P P 2 ^ 


- ? (-1)” -^AA2^tanh(u8j,t) + AA4^tanh(v6^t)] 
n*l 


,n ®n 


aa 


Case 2 . Roots Real and Equal in Pairs (±u, v > 0) 
Syimetric x, symmetric y (SS); 


SSr + (E1A„) ? (E1B„„)SS2„ » (E1A„) {(E1C„) + ? (tlD„„)} 
P P p.] ' pn' n ' p' ' p' pn' 


(FIAJ ? (F1B„„)SS1 * SS2„ - (FlAl {(FlC ) » ? (F1D„J) 


n»l 


pn' n 


n=l 


pn- 
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where 


ElAp = 



r 2 

ua„£tanh uci„Jl 

2a.^t 

p 

1 _,P P 

tanh(ua + ua «. 


ElBp^ = 4{-l)P^’' 


( ue ,)" 

-— ? ■ tanh(u6„t) 


Ell 


£1 


[tanh(uB^t) + uS^tJ 


E'S * ■ 


a„Jltanh(ua„j!,) 

- P P c 

'l+2p ’ ttanh{uapS,) + ucpil] 5-6p 


E1D -3(1)P'*'^ ^ . 

C(uSJ^ + a^] u6^Jl[tanh(u6 t) + u6, t] 


pn 


n ' ri 


Ell 


\ [( uS ^)^ + a :]/ 


pn 




tanh(uB^t) + uS^t 


FUp » 


1 


Z&ll 


1 - 


uSpt tanh^uBp" 
tanh{uBpt) + u 6 pt J 


F1B.„ = 4(-l)P‘'" 5 ^ tanh(ua 0 

P" [(uc/ + b 2 ] " 


— -1 
tanh(ua„Jl) + ua„il 
n n 


ncp > - 


Bpt tanh(uBpt) 

^3+4p ■ [tanh(uB-t) + u6_t] ^7-8p 


P ' P 
^3 


F 11 


fid ■ 2 (- I )'’*" j = S . 

"" [(uc„)2 + S‘] ud^[tanh(ua„i,) ♦ uc.I] 


n 


n ' n 
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/ [(uoj^ - S?] \ 

FIT _ » (2 ^ ) tanh(ua £) + ua J, 

"" V [(ua„)2 * tl: } 


Also, 


“ tanh(uap£) + uOpJl 

“ tanh(u8pt) + u8pt 
''xss ° It '^1+20 


S5.6P - S51p t»nh(uapl) 


— ~T' S7 On - SS 2 _ tanh(u8nt) 
2UB^H P P 


i'6 

- ? (“1)*^ -r^ [SS2„ tanh(u8-t) 

n»l ^ " n 

+ SS4^ {tanh(u8^t) + u6^t}] 

•^yss “ ir ^3+4o 

- ? (-1)" -r^ [SSI tanh(ua £) 

n=1 ^ " n 

+ SS3^ {tanh(ua^Jl) + ua^ <1}] . 


S>'mmetric x, antisymmetric y (SA): 


S*’p + (E2Ap) (E2Bp„)SA2^ - (E2Ap) {(£20,,) * (E2Dp„)} 

(F2Ap) (F2Bp„)SAl„ ♦ SAE^ • (FEA^) {(F2C^) * (F20p„)) 
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where 


E2Ap 




1 

* 

Uttpi tanh (uOpA) 
Ltanh(uapJl) + uapJlJ 


E2B„„ = 4 E21„„ tanh(u6„t) 
pn pn ' n ' 


ttpjl tanh(uap£) 

^ [tanh(ua i) + ua i] ^5-6p ' ^l+2p 


* 2 
pn 


[(u 6^)2 + a^] P" 


'3-4n 


P" [(uej^ * "®n‘ \r(u6j^ + c?] / 


[(U6„)‘ * o.p 




tanh(u6^t) 


' ^^p 2Sp£{u6pt - [1 + u6pt tanh(u3pt)]tanh(u6pt)} 


ua^S 

F2B^_ = 4(-l)P'''’ p-?- tanh(ua £) 

P" lim/ + Bp " 


F2Tn„ 

1 _ pn 

[tanh(ua^^«,) +' ucx^2-] 


‘■'^^p " ^7+8p ’ Bpt ''5-60 


(-l)P . . tanh(uB„t)] 

+ _ . Cr ■*■ 


3„t 4-4p 




F21 


6n 
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F21 


pn 


tanh(ua^i,) - 


2 ( ua ^)2 

5 s- tanh(ua„Jl) + ua„£ . 

C (%)^ - 8 ^] 


Also, 


SA3 


■1 


p ~ [tanh(uapi) + 


SAl„tanh(ua 2) + 

P P Zua^t 


SA4 


-1 


p U6pt 


SA2„tanh(u6„t) + 

p P 


2 b ;£ 


5 - 3 - 4 o 

7 sa 4 £t 


Antisymmetric x, symmetric y (AS): 


AS! + (E3AJ ? (E3B_JAS2„ = (E3AJ {(E3Cj + ? (E3D_J) 

K P pn n p p ^_1 pn 

( F 3 Ap ) ( F 3 Bp „) ASl „ + AS 2 p = ( FSA ^) {( F 3 Cp ) + F 30 p „)} 

where 


p, . ^ 

^''^p 2apt{uap2. - [1 + UttpJl tanh(uapJ2,)]tanh{uapJi) } 


2 

„ . ^ u 8„a_ 

P ^ tanh(uB^t) 


E3B =4(-l) , ^ 

[( U 8 J " t o,b 


E 31 n 

pn 

[tanh(uBpt) + u8^t] 


f.i^P 1 + uot.i tanh(ua_e) 

P 3 C = S 4 * — C 4 * — P P r 

“S '5*6p apt <^7-80 apt ^l-2p 


1 ^ 


1 



T 
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. a E31 

l[(u6„)^ ♦ 0.2] UanMu8„t) ♦ uB„tJ 

2(u6n)^ 

* tanh(u8-t) 5 x- tanh(u6„t) + u3„t 

I’" " [(u8„)2 ♦ C.2] 


F3A, 


1 


2Spi 


1 - 


u0pt tanh‘^(u6pt) 

[tanh(u6„t) + u0„t] 


= 4 F31_„ tanh(ua„£) 
pn pn ' n ' 


8pt tanh(uBpt) 

''^^p “ Ltanh(u6pt) + uBptJ ^7-8p ‘ ^3+4p 


F30 


pn 




(■1)P*" , F3, 

[{ua„)2 * 8p 


pn 


'l-2n 


F31p„ = 2(-l) 


p+n 


%2p 


i[(uap)^ + lC(uaJ^ + Sp] 


(uaj- 


-1 


Also, 


AS3. 




ASl tanll(ua„l) ♦ 


AS4_ = 


-1 

tanh{u8pt) + uSpt 


AS2„ tanh(ua^t) + 

P P 2u6?£ 


, = 1-20 
>as 4it ■ 


^7-8n 


tanh(ua i) 
n 
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Antisyrmetric x, anti synmetri c y (AA): 


AAl + (E4AJ ? (E4B„JAA2„ - (E4AJ {(E4Cj + ? (E4D„J} 


n«l 


pn' n 


n=l 


pn' 


(F4AJ ? (F4B„JAA1 + AA2„ = (F4Aj {(F4C ) + ? {F4D„J> 


n=1 


pn' n 


n*l 


pn' 


where 


E4A, 




2 ua;t 


1 - 


(1 + uajl tanh(ua_«.)) 


.. — tanh(ua„t) 

ua_£ p 


(u6„)^ E41 „ 

FdR - 11 5 7 ~ tanh(u6„t) 1 — 

E’Bpn - [(ue„)2 + o^] " " 


1 + ua fi, tanh(ua il) 

2 ^ E_ s - Cr 

p QipJl l-2p 5+6p 


E4D 


(ue„)‘ 


E41 


£n 


pn “ [(uB„? + Qip] u6^n 


E41 „ = u6„t + 
pn n 


- 4: - 


tanh(ua^t) 


F4A, 


1 




[1 + u6 t tanh(u6 t)l 
Wt ^ t»r,h{p 6 pt) 


2 

( ua_ ) a_ 

F4B = 4(-l)P'^" 5 V tanh(ua 1 ) 

* b‘] " 


F41 


1 - 


pn 
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1 + uB-t tanh(u6_t) 

E ! — 2 s 

uBpt ^3-4p 


^7+8p 


F4D = 2(-l)P'^ 2 S, 2n 


F41 „ * ua J2, + 
pn n 


1 - 


[(ue.„)^ - B?] '' 

[(“«„)'* 6^] J 


tanh(ua^£) 


AU?. 


AA3_ = 


P Uttp^ 


AA1_ tanh(ua„Jl) + 

" " 2o?t 

P 


AA4p 


uBpt 




^ 5 - 1^0 

~JT~ 


? 

n=l 


tanh(uB-t) + 

(-1)” ^ [AA2^ tanh(u6„t) + AA4^ uB^t] . 
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APPENDIX D 

INCLUSION OF DISPLACEMENT BOUNDARY CONDITIONS IN THE ANALYSIS 

To include displacement boundary conditions in the exact solution 
analysis, the problem is reformulated using the displacements u (for 
the x-direction) and v (for the y-direction) as the variables of interest. 
All of the assertions made earlier as to the character of the problem 
(plane stress, rectangular slab, etc.) are still considered to hold 
true. 

At this point, the normal and shear strains are defined by the 
relations 

S ' “’X • s “ ’'y • V ' “’y * "’X • 

Also, the stiffness version of Hooke's Law is given as 
^x = ^llS ^12^ ^ ^16% 

^y " ^12^x ^22^y ^26"^xy (D-2) 

"xy ^ ^16^x ^ ^26^y ^66^xy ’ 

where the stiffnesses Q.. are as defined in Jones [10]. Specializ''’i 

* J 

the problem to orthotropic materials having their principal axes ali^ined 
with the x- and y- directions of the slab makes 

‘^16 ' ^26 " ° • 


The equilibrium equations then become 
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^’xx ^66 ^’yy ^^12 ^66^ '^’xy “ ° 

•^ee ''’xx ■*■ ^22 ''’yy ^*^12 ^66^ “*xy * ° ' 

Differentiating equation 0-3 with respect to x gives 


(D-3) 

(D-4) 


0, , u, + Q/rr u, + (Q-i/i + Qt/;) Vj = 0 . 
’XXX ^66 ’xyy '^12 % 6 ‘ ’xxy 


Differentiating equation 0-4 with respect to y gives 


^66 ''’xxy ^22 ''•yyy ''' ^^12 ^66^ ^’xyy * ° 


(D-5) 


or 


''•xxy “ ■ Q.. ''’yyy 


Q22 (^12 ^66^ 


^66 %6 
Combining equations D-5 and D-6 gives 




(D-6) 


‘'’xxx * ^66 ‘'’xyy 


(^12 ^66^ 
^66 


(D-7) 


^22 ''’yyy ^^12 ^ ^66^ ‘'’xyy^ " °’ 


Differentiating equation D-3 twice with respect to y and rearranging 
gives 


'n 
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''’xyyy “ ’ (Q^2 '^ ‘^66^ ^’xxyy " (0^2 ^ ^66^ ^’yyyy 
This can be combined with equation D-7 and manipulated to give 


(D-8) 


u, +Au, +Bu, =0, 
’xxxx xxyy yyyy 


(D-9) 


where 


Wat 

, - - i 



r 
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= *^22 _ ^ 12 ^ ^^12 ^ p ,J_ _ ^^ 12 , 

^ ^ 11^66 ^11 ^ *"12 


B 


= ^22 , ^2 
^ ^1 


Similar manipulations also yield 


V, + 

xxxx 


^ '^’xxyy ® '^’yyyy ~ ^ 


(D-10) 


with A and B defined as above. 

At this point, the equations D-9 and D-10 have the same form as 
equation 4, with u and v, respectively, in place of 4). 

With suitably complete solution forms for u and v, equations D-3 
and D-4 can be used to relate the coefficients of one of the solution 
forms to the coefficients of the other. This done, the solution 
becomes a problem in u or v alone. Next, any stress boundary conditions 
must be expressed in terms of displacements. Following the same basic 
path as the development for stress boundary conditions then leads to 
the displacement form of the solution. 




i 



